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Abstract 

We introduce the concept of multiplication matrices for ideals of pro- 
jective dimension zero. We discuss various applications and in particular, 
we give a new algorithm to compute the variety of an ideal of projective 
dimension zero. 

1 Introduction 

Eigenvalue methods to compute the variety of an affine zero-dimensional ideal 
has become an active area of research EJ \W[ \T7\. Recall that an ideal is 
defined to be of dimension zero when the corresponding variety is finite. The 
eigenvalue methods use both algebraic and numerical algorithms and the key 
is a nice one-to-one correspondence between the points on the variety and the 
eigenvectors to so called multiplication matrices. 

The notion of zero-dimensionality has a meaning also in projective space. 
Over an infinite field, an ideal of projective dimension zero is an ideal whose 
variety consists of a finite number of projective points. Below we will give a 
general definition, valid also for finite fields. We show that it is possible to 
define multiplication matrices with respect to ideals of projective dimension 
zero. Our main result is that the one-to-one correspondence mentioned above 
also holds in the projective setting, giving a new method to compute the variety 
of an ideal of projective dimension zero. 

In order to define the projective multiplication matrices, we need to choose 
appropriate vector space bases for the graded pieces of the quotient ring S/I, 
where S is a polynomial ring over some field and / is an ideal of projective 
dimension zero. 

The usual choice of a vector space basis is the set of residues to the comple- 
ment of the initial ideal of I (with respect to some monomial order) . Our choice 
of bases differs from the usual ones — in general we consider non-monomial 
Ik-bases. This choice of bases happen to give a fast normal form algorithm for 
high degree elements, which outperforms the usual reduction method based on 
Grobner bases. 

Moreover, we discuss how our approach can be used to compute vanish- 
ing ideals of projective points. We give an alternative version of the graded 
Buchberger-Moller algorithm and give a fast algorithm for computing projec- 
tive separators. 



2 Notation and preliminaries 



Throughout the paper, let k be a field and let S = k[xo, ■ ■ ■ , x n ] denote the poly- 
nomial ring in n + 1 variables. Recall that the Hilbert series of a graded ring 
R = Ro © Ri © Ri © • • • is the power series Hs(i?, t) — dimk(-Ro) + drmjt(-Ri)f + 
dimt(R 2 )t 2 + ■ • ■ . An ideal / is of projective dimension zero exactly when 
R = S/I is graded and satisfies dimj^i?.;) = m for some m > and for all 
i sufficiently large. The least i such that dim^iii) = dimk(-Rj+i) = • • • is called 
the postulation number and is denoted by post(i?). We say that R postulates 
in degree post(iZ). 

When / is an ideal of projective dimension zero and R — S/I, we say that 
R is a ring of projective dimension zero. For an element a of S, we write [a] to 
denote the equivalence class in R containing a. 

By V(I), we denote the variety of / with respect to the algebraic closure 
k of k, so |V(/)| < The number of projective points in V(I) counting 

multiplicity equals dim k (i? post ( fl) ). 

Let / = Q\ D • • • D Q s be a minimal primary decomposition of /. Let Pj = 



s/Qi- When no P\ equals the unique graded maximal ideal m = (xq, . . . ,x n ) 

of S, we say that / is unmixed. When / is unmixed, the Hilbert series of R is 
strictly increasing until it reaches degree post(i?). 

When / is mixed, we will write the primary decomposition as / = Qi n 
■ • • ("I Q s H Qs+i and assume that m = -J Q s +i ■ The Hilbert series of R when 
/ is mixed does not behave nice in general (it can have an arbitrary number of 
valleys for instance). 

We need a notation to drop the primary component Q s +i in the case when 
/ is mixed, so we define I u = Q±H - ■ -r\Q s . This means that when / is unmixed, 
we have the identity I = I u . We extend this definition so that R u = S/I u . 

Example 2.1. Let I = x\x\ — x\, x\x\ — x\. Then I = (x\ — X2)n(a:|)n(a;f, x%) 
Thus, Qi = (xi - x 2 ), Q 2 = {x\) and Q 3 = (xf ,x\). We have Pi = y/Q~l = Q u 
so V{Px) = {(1 : 1)} and P 2 = VQl = (x 2 ) so V(P 2 ) = {(1 : 0)}. Finally 
m = V^3- Thus, I u = (xi - x 2 ) n {xl) = {x y x\ - x%). We have Rs{R u ,t) = 

1 + It + 'it 2 + M 3 H , while Hs(i?, t) = 1 + 2t + 'it 2 + At 3 + 4< 4 + it 5 + it 6 + ■ ■ ■ , 

so post(i? u ) = 2 and post (R) = 5 (It does not in general hold though that 
post(i? u ) < post(i?).) The point (1 : 1) has multiplicity one, while the point 
(1 : 0) has multiplicity two (it follows from the fact that (x 2 ) has a chain of 
primary ideals of length two) . 

We have a one-to-one correspondence of prime ideals Pj generated in degree 
one and points on V(I). If V(I) = V(I), it follows that \V(I)\ = s and that 
Pi, . . . ,P S are all generated in degree one. 

The concept of non-zero divisors is of particular importance in this paper. 
Recall that I is a non-zero divisor on the S'-module M if Im = implies that 
1^0. In Example 12.11 [xi] is a non-zero divisor on R u , while [x 2 ] is not. The 
existence of non-zero divisors is connected to Cohcn-Macaulayness and also to 
the primary decomposition in the following sense for a ring R = S/ 1 of projective 
dimension zero. 

/ is unmixed R is Cohen-Macaulay R contains a non-zero divisor 

Although R lacks non-zero divisors when / is mixed, we will show that there 
exists a minimal degree d so that R% © Rd+i ® " ' anc ^ © Rd+i © • • • are 
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isomorphic as S'-modules. This degree d will equal max(post(i?), post (.ft")). To 
simplify notation, we will denote max(post(ft), post(ft u )) by nz(ft). When the 
context is clear, we will omit ft and only write nz. It follows that every non-zero 
divisor on R u is also a non-zero divisor on ft nz © R n z+i © • ■ • . In Example 12. 1[ 
nz(ft) = 5 and [x{\ is a non-zero divisor on R§ © Rq • • • . 

Suppose that V and W are two k-spaces of dimension m and m' respectively. 
Let ei, . . . , e m and /i, . . . , / m ' be k-bases of V and W, respectively. Let ^ be a 
k-linear map from V to W and let be the m x m' matrix, whose i'th row is 
the coefficient vector (ci, . . . , c m >) and where 0(e,) = ci/i + ■ • ■ + c m 'f m '. Notice 
that is the transpose of the standard matrix representation of <f>. 

When the map <f> is defined on an a finite dimensional algebra by v i— > fv 
for an element / in the algebra, then the matrix is called the multiplication 
matrix with respect to /. We also denote A^ by Af 

Example 2.2. Let I = (x% — 1, x| — X2) C Q[xi, xq\. The ring Q[xi, X2]/I is of 
affine dimension zero and a vector space basis for the quotient ring can be chosen 
as [l],[x 2 ]. In this algebra we have [xi][l] = [1], [x±] [x 2 ] = [x 2 ], [x 2 ] [1] = [^2] 
and [x2\ [xz] = [X2], so the multiplication matrices with respect to [xi] and [x 2 ] 
equals 



When it comes to rings of projective dimension zero, we have a problem 
since these rings are infinite dimensional. In order to overcome this problem, 
we use the graded property among these rings. So let <ft be defined from Rd to 
■Rd+l/i by multiplication by a form / with respect to the bases e\, . . . , e m and 
fx, ... , f m i for Rd and respectively. Then Af (or A$) is the projective 

multiplication matrix in degree d with respect to /. We will later show that it is 
possible to choose bases such that the projective multiplication matrices agrees 
for all degrees greater than or equal to nz(_R). 

Example 2.3. For the ring R u from Example 1 2. 11 we can choose [xf\, [xix 2 ], [x 2 ] 
as a k-basis in degree two and [xf\, [x\x 2 ], [x 2 ] as a k-basis in degree three. If 
<j) denotes the map from i? 2 to R% induced by multiplication by [x 2 ], then the 
projective multiplication matrix of degree two with respect to [x 2 ] equals 



since [x2][xf] = 0[xf] + ljxfa^] + 0[x^], [x2][xiX2] = 0[xf] + 0[a^x 2 ] + l[x|] and 



3 Projective multiplication matrices 

We will use the fact that rings of projective dimension zero postulates in order 
to define the projective multiplication matrices. 




It is customary to write A Xi instead of Ar^i and we will do so in the sequel. 




[x 2 ][xf] =0[x?]+0[x?x 2 ]+l[x!]. 
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3.1 Non-zero divisors of degree one 

If pi = (pio : ■ ■ • : pi„) is a projective point with respect to the coordinates 
xq, . . . ,% n , then we cannot define evaluation on a form in a unique way since 
Xpi — pi for non-zero A. For practical purposes it is however good to define 
the evaluation in a unique way and we will do this by simply fixing a repre- 
sentation of each projective point. So we associate to each projective point 
Pi = {Pio Pin) an affine point pi = (pio, ■ ■ ■ ,Pin) and now we can define 

the unique evaluation as Xq° ■ ■ ■ x"™ (pi) — Xq° • • • x"™ (pf) = p" ° • • • p^ ■ This 
way of thinking of evaluation is implicit in [18] . 

With this evaluation method it follows that if h and h are elements of S 
and [li] = [I2] in P, then h(p) = hip) for all points p on V(I). This property 
makes it possible to define evaluation on elements in R by [l](p) = lip), where 

P eV(l). 

We will need a notation for evaluating a set of points on a set of ele- 
ments. If P = {pi, . . . ,p m } is a set of projective points, we write f(P) = 
(/(pi), • ■ • , f(pm))- If F = . . . , /,} is a set of forms in S, then F(P) is 
defined to be the (s x m)-matrix whose z'th row is fi(P). 

Lemma 3.1. Let I be an ideal of projective dimension zero. Suppose that 
V(J) = V(J). That [I] is a non-zero divisor on R u is equivalent to l{p) ^ for 
allp£V(I). 

Proof. Suppose that lip) ^ for all p G V(I). If [a] G R u is such that [a] ■ [I] = 
in R u , then a • I G I u , so that (a • l){p) = a(p) • lip) = for all p e V(I). Hence 
a(p) = for all p € V(I). Thus a £ I, so [I] is a non-zero divisor on R u . 

Suppose instead that [I] is a non-zero divisor on R u . Let q be an arbitrary 
point in V(I). Let [Q] be an element in R u such that <5(?) 7^ and Q(p) = for 
p G V(7) \ {g}. (The element Q is called a separator for the point q with respect 
to V{I). Separators exists, see for instance pQ or Section 14.4.11 ) Suppose that 
l(q) = 0. Then Qip) -lip) = for all p G V(I), so that [Q]-[l] = 0. Since [Q] ^ 0, 
it follows that [I] is a zero-divisor, which is a contradiction. Hence l(p) ^ for 
aJlpGV(J). □ 

Proposition 3.2. Suppose that k contains at least |V(/)| elements. Then R 
has a linear non-zero divisor if and only if I is unmixed. The requirement on Ik 
is sharp in the sense that ifh contains |V(7)| — 1 elements, then there exists an 
unmixed ideal I such that R lacks linear non-zero divisors. 

Proof. If I is mixed, then R does not contain non-zero divisors. So suppose 
that I is unmixed. Let s = \V(P)\ and let I = Q\ n • • • PI Q s be a primary 
decomposition with respect to S. Let Pi = \fQi. Each Pi is generated in 
degree one. The set of zero-divisors in R equals the union of the residues of the 
Pi'.s (Proposition 4.7 in 2J). So if we let NZ = Sx \ UP,, then the set of linear 
non-zero divisors in R is the residues of NZ. Let NZ — Si n NZ. Then the 
linear non-zero divisors of R is the residues of NZ. Let r\i — Pj n S\. 

Suppose that k is infinite or finite and contains |k| > s elements. Suppose 
that Vi G Si but vt £ r\i U • • • U rji for i < s (clearly Vi exists). If Vi (£ rjj+i, then 
let Vi+i — Vi. Otherwise, take an element Wi G rji such that Wi ^ r?i+i (such 
an element must exist, since we assume that Vi ^ f?i+i)- The element Vi + aw, 
does neither belong to r\i nor to r\i+i for any non-zero a G k. Pick oti G k\ {0}. 
If Vi + ctiWi G Tji—i, then Vi + awi (fc rji—i for all a G k \ {ai}, since otherwise 



4 



we would have a(vi + a±Wi) — a%(vi + awi) = (a — ai)vi G ift— 1, which is a 
contradiction. So pick a 2 G k \ {a{\. Clearly rji + a 2 wi ^ rji-i U rji U ^i+i- It 
is clear that we can continue in this way provided that there is at most i — 1 
non-zero elements in k. Since i ranges from 1 to s, this construction uses at 
most s — 1 non-zero elements in k. 

Suppose instead that k is a field with elements {ag, a\, a 2 , . . . , a s „i} with 
«o = and a± = 1. Consider the points p% = (1 : : • • • : 0),P2 = (1 : 1 : 
: ••• : 0),p 3 = (1 : a 2 : : ••• : 0),...,p 3 = (1 : a s -i : : ■ ■ • : 0) and 
Ps+i = (0 : 1 : : • • • : 0) in P n (k). Let /(pi) be the vanishing ideal with 
respect to pi, which is prime. Let / = C\iI(jPi) and let I = bgxo + • • • + b n x n 
be an arbitrary linear form. If I is non-zero on pi, . . . ,p s , then bo + b\ao 7^ 0, 
bo + bidx 0, . . . , bo + &ia s _i ^ 0. Now b\ must equal zero, since otherwise we 
would have bo+b\a,i ^ bo+bidj when i ^ j and thus, by the pigeonhole principle, 
bo + bidi = for some i, which contradicts the assumption that l{pi) ^ 0. But 
if b\ = 0, then l(p s+ i) = 0. Thus, the ring lacks linear non-zero divisors by 
Lemma 13.11 

□ 

Example 3.1. Let I — (xo +Xi,xq +X2) H (xo +xi,xq + 2x 2 ) fl (xq + 2x\,xq + 
2x 2 ) C Z 3 [xo,aJi,aJa]. We have V(i) = {(1 : 2 : 2), (1 : 2 : 1), (1 : 1 : 1)}. To 
compute a non-zero divisor, we start by computing v\. We pick an element of 
degree one in (xq + Xx,Xo + X2), say xq + x\. By changing one of the coefficients, 
we can assure that this element is not in (xo+x\, 2:0+2:2), so we let v\ = xq+2x\. 
Since V\(p2) ^ we let v 2 = v\. But ^2(^3) = 0, so we look for an element w 2 
which is in (xq +xi,xq + 2x 2 ) but not in (xo + 2xi, xq + 2x 2 ). It is clear that we 
can find such an element by going through the generators of (xq +x±,xq + 2x 2 ) 
until we find an element which is not in ^3. Indeed, w 2 = xq + x\ is such an 
element. We have v 2 +w 2 = 2xq. Since 2xo(pi) 7^ we can use [2xo] (or rather 
[2:0]) as a non-zero divisor. 

When k is finite, it is an interesting question to determine, given a degree d, 
the maximal number of points allowed to guarantee the existence of a non-zero 
divisor of degree d. Using a result due to Chevalley, one can show that there is 
always a non-zero divisor of degree n, even if all points in P n (k) are considered. 
However, we will only use non-zero divisors of degree one in this paper and this 
problem will be dealt with in a separate paper [13] . 

Proposition 3.3. Let I be an ideal of projective dimension zero. Suppose that 
k contains at least |V(/)| elements. Then there exists a linear form I G Si such 
that £ : Rd — > Rd+i, [a] | — > [l][o] is onto, for all d > nz(i?). 

Proof. If R = R u , then, by Proposition 13.21 R contains a nonzero-divisor of 
degree one which has the desired property. Otherwise, the maximal ideal is 
associated to /. Thus, the primary decomposition of / can be written as I — 
J D Q, with y/Q = m. Let d > nz(R). Then dimi s (i? ( j) — dim^i?^), which is 
equivalent to dimj^Jd) = dim^ J n Q)d- Since (J n Q)d C Jd 7 it follows that 
Jd = (J H Q)d- So that I is a nonzero-divisor on ®i>dR^d implies that / is a 
nonzero-divisor on (Bi>dRd- d 

Suppose that {[ei], . . . , [e m ]} is a basis for the k-space Rd, d > nz(i?). Since 
the map induced by / is onto, the set {[eiZ], . . . , [e m Z]} forms a k-basis for Rd+i- 
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In general, for any positive integer i, {[e\l l ], . . . , [em? 1 ]} can be chosen as Ik-basis 
for Rd+i- This shows the following important theorem. 

Theorem 3.4. Let I be an ideal of projective dimension zero. Let {ei, . . . , e m } 
be a k-basis for R nz . Suppose that R u has a non-zero divisor [I] of degree one. 
Then {[eiT], . . . , [e m P]} is a k-basis for i? n z+i, for all positive integers i. The 
matrix multiplication matrix Aj representing the map R n z+i — > Rnz+i+i, [a] >— ► 
[xjd] with respect to the bases above, is independent of the choice of i. 

3.2 An affine connection 

To a ring R of projective dimension zero, we will now associate an affine ring of 
dimension zero — i?*, whose multiplication matrices coincide with the projective 
multiplication matrices of the projective ring. In fact, the zero-dimensional ring 
i?* is simply R/([l] — [1]). The key is the following lemma. 

Lemma 3.5. Suppose thatk contains at least \V(I)\ elements. Then there is a 
linear change of coordinates T and a variable Xi such that T(xi)(j>) ^ for all 
points p £V(I). 

Proof. Let I be the form from proposition l3.3l We can write I = boXo + - 1 ■+b n x n . 
Some coefficient is non-zero, say 6, 7^ 0. Let T(xi) = I and let T(xj) = Xj if 
j^i. □ 

Remark 1. With T as above, if we let yo — T(xi), y\ — T(xq), 1/2 — T(xi), . . . , 
yi = T(xi-i),y i+1 = T(x i+ i), . . . ,y„ = T(x n ), then each point on V(I) with 
respect to yo, . . . ,y n can be written as (1 : a% : ■ ■ ■ : a n ). 

Let J be an ideal of affine dimension zero and let J = q±n- ■ -ng s be a minimal 
primary decomposition. The multiplicity of a point p S V(J), belonging to the 
primary component is defined as the length of qi. 

The same definition holds for projective points, that is, if I is an ideal of 
projective dimension zero and I = q\ PI • • • PI q s , then the multiplicity of a point 
p on V(I) belonging to the primary component y/qi, is defined as the length of 
Qi- 

These two multiplicity definitions are connected in the sense that if p = (1 : 
a\ : ■ ■ ■ : a n ) is a projective point with multiplicity r in V(I), then (ax, . . . , a n ) 
is an affine point with multiplicity r in V(I + (yo — 1)). This is a standard result 
and treated in [TH] and [H] for instance. 

Lemma 3.6. Let I be an ideal of projective dimension zero. Suppose that there 
exists an I £ S\ such that [I] is a non-zero divisor on R u . Put i?» = R/([l] — [1]). 
Let {[ei], . . . , [e m ]} be a k-basis for Rd, for d > nz(i?). Then {[[ei]], . . . , [[e m ]]} 
is a k-basis for R*, where [[ ]] denotes an equivalence class in R mod [I] — [1]. 

Proof. Since the k-dimension of R/([l] — [1]) and i? nz is determined by the sum of 
the points counting multiplicity, we have dimk(i?/([Z] — [1])) = dimj^-Rnz) = m. 
Since [I] is a non-zero divisor on the basis i? nz , the vectors [[ei]], . . . , [[e m ]] are 
linearly independent. □ 

Proposition 3.7. Let I be an ideal of projective dimension zero. Suppose thatk 
contains at least \V(I)\ elements. Let m = dimk(_R nz ). Put i?* = R/([l] — [1]) = 
S/(I + ([I] — [1])), where I is a non-zero divisor of degree one on R u . Let 
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{[eiP], . . . , [e m Z*]} be a k-basis for Rd+i and let {[[ei]], . . . , [[e m ]]} be a k-basis 
for i?*, where [[ ]] denotes an equivalence class in R mod [I — 1], Then there 
is a change of coordinates such that the multiplication matrices with respect to 
Xi, . . . ,x n coincides for i?„ and R. 

Proof. By a change of coordinates in accordance with Remark[T] we may assume 
that I = y,Q. Let Ak = be the projective multiplication matrix of R with 
respect to yk, such that [y^jfej = aii[yo e i] + • • • + o,i m [yoe m ]. It follows that 

[[2/fe]][N] = o»i[[ei]] + hOi m [[e m ]]. By Lemma EH1 {[[ei]], . . . , [[e TO ]]} is a 

k-basis for Hence Ak is the multiplication matrix of i?* with respect to yk- 

□ 

3.3 Computing the variety from the projective multipli- 
cation matrices 

We now state the result of Moller and Stetter in the afhne setting. See [TH] for 
a proof. 

Theorem 3.8. Let I be a zero dimensional ideal. Let {[ei], . . . , [e m ]} be a 
k-basis of k[x\, . . . , x n ]/I. Let A±, . . . , A n be the multiplication matrices with 
respect to this basis. Let r — |V(/)|. Then there are exactly r common (right) 
eigenvectors for the matrices A\, . . . , A n and they are (ei(pi), . . . , e m (pi)) t for 
i = 1, . . . , r. Let Ay denote the eigenvalue of Aj corresponding to the eigenvector 
{ei{pi),...,e m {p i )) t . Then pi = . . . , X im ). 

We have an almost identical theorem in the projective setting. 

Theorem 3.9. Let I C k.[xo, ■ ■ ■ ,x n ] be an ideal of projective dimension zero. 
Suppose that k. contains at least |V(/)| elements. Let {[ei], . . . , [e m ]} be a k- 
basis of R nz and let I be a linear form such that {[eil], . . . , [e m Z]} is a k-basis 
for i?nz+i- Let Ao,...,A n be the projective multiplication matrices with re- 
spect to this basis. Let r — |V(7)|. Then there are exactly r common (right) 
eigenvectors for the matrices Aq, . . . , A n and they are (ei(pj), . . . , e m (pj))* for 
i = 1, . . . , r. Let Ay denote the eigenvalue of Aj corresponding to the eigenvector 
(ei(pi), . . . ,e m (pi))*. Then p t = (X i0 : X a : ■ ■ ■ : X im ). 

Proof. By Proposition 13. 3[ there exists a linear form I = bgXQ + • ■ ■ + b n x n 
such that {[eiZ], .. . , [e m Z]} forms a k-basis for R nz+ i. By Remark [l] there is a 
change of coordinates such that the multiplication matrix Bi with respect to yi 
satisfies B = b A + • • • b n A n and B\ — A ,B 2 = A\, . . . , Bi = Aj_i, B i+ i = 
Ai+i, . . . , B n = A n . 

By Proposition 13.71 the projective multiplication matrices B\, . . . , B m of R 
with respect to y\, . . . , y m agree with the multiplication matrices for R/(l — 1) 
with respect to y\ , . . . , y m . 

Since the multiplication matrix with respect to yo is the identity, a common 
eigenvector for B\ , ... , B n is also a common eigenvector for Bo, ... , B n and vice 
versa. But by linearity, v is a common eigenvector to Bq, . . . , B n if and only if 
v is an eigenvector to Aq, . . . , A n . Hence, the set of common eigenvectors for 
A ,. . .,A n equals ei(pj), . . . ,e m (pi), for i = 1, . . . , r, by TheoremOTEl □ 

To determine the multiplicity of a point p e V(I), one can use the result 
of Corless et al in [5]. The method goes as follows. Let A be a generic linear 
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combination of the multiplication matrices. Let A be the eigenvalue of A with 
respect to e(p) (clearly e(p) is an eigenvector of A). Then the multiplicity of p 
equals the algebraic multiplicity of A. There are also direct methods which one 
could use, see for instance [E] and [T7] . 

Example 3.2. The elements 

/i = xz + yz - z 2 

f2 = x 2 -V 2 + 2yz - z 2 

/3 = xy - y 2 + yz 

generates an unmixed ideal / of projective dimension zero in <C[x, y, z\. Choosing 
[x] , [y] , [z] and [y 2 ] , [yz] , [z 2 ] as bases in degree 1 and 2 respectively, we see that 
neither [x], [y] nor [z] serve as nonzero-divisors. Indeed, if we let M x ,M y and 
M z denote the multiplication matrices from i?i to i?2 with respect to the bases 
chosen above, we compute 



M x = 




M z = 




and we can see that all the matrices have a nontrivial kernel. However, M y + M z 
has full rank which is equivalent to [y + z] being a nonzero-divisor. Hence, if we 
use + [y{y + z)], [z(y + z))} as a k-basis in degree two, we can construct 
the projective multiplication matrices A x ,A y and A z . From these matrices the 
solutions can be read off. Now [x(y + z)] = [y 2 ] — 2[zy] + [z 2 ], [y(y + z)] — 
[y 2 ] + [zy], [z(y + z)] = [yz] + [z 2 ] by making use of the multiplication matrices 
above. Thus, with 



T 




we have 

A x = M x (T t ) 



and similarly for A y and A z , so that 



A x — 



and 




Common eigenvectors for the matrices are (1, 1, 0), (1, 0, 1) and (0, 1, 1). The 
eigenvalues corresponding to (1, 1,0) are 1, 1,0 for A x ,A y and A z respectively. 
Likewise, the eigenvalues corresponding to (1, 0, 1) are 1, 0, 1 and the eigenvalues 
corresponding to (0, 1, 1) are 0, 1, 1. Thus, V(I) = {(1 : 1 : 0), (1 : : 1), (0:1: 
!)}• 

Notice that since nz(i?) = 1, we can also use the correspondence between 
eigenvectors and the k-basis to obtain the points. Indeed ([x](pi), [y](pi), [z](pi)) 
— (1, 1,0), thus we have pi = (1 : 1 : 0), etc. 
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4 Applications and computational aspects 



A convenient way to think of a ring R — S/I of projective dimension zero is as 

R = R © Rl © • ■ ■ © Rnz 

together with the linear map I and the multiplication matrices A%, . . . , A n . We 
write this information as a triplet (R,A,l). The Ik-dimension of the graded 
pieces of R describes the configuration of the points and also tells whether or 
not the maximal ideal is associated, while the multiplication matrices encode 
the variety as a set. In Section 14.11 we will see that we obtain a fast normal 
form algorithm by using the triplet. With this perspective, the classical way of 
determining a Grobner basis for / misses a lot of information about the ring. 
It also turns out that we compute unnecessary data. For instance, a Grobner 
basis for the ideal I = (xz + yz — z 2 , x 2 — y 2 + 2yz — z 2 ,xy — y 2 + yz) from 
Example 13.21 with respect to x > y > z and DegRevLex is (xz + yz — z 2 ,x 2 — 
y 2 + 2yz — z 2 ,xy — y 2 + yz,y 2 z — yz 2 ). Since nz(i?) = 1, we only need to consider 
the Ik-spaces R\ and i?2 to determine the variety, and for this purpose, the term 
y 2 z — yz 2 in the Grobner basis is superfluous. In Section l4~2l we will show that 
the maximal degree of a term in a Grobner basis is max(nz(i?), m). Since it is 
enough to compute up to degree nz(i?) in order to determine the variety, this 
indicates that Grobner techniques are not always optimal. Unfortunately, it is 
hard to detect nz(R). 

4.1 Computing normal forms with respect to (R,A,l) 

As an application of the multiplication matrices, we obtain a fast normal form 
algorithm for high degree elements of S. Suppose that we have a normal form 
algorithm Nf(*, B) for elements of degree less than or equal to nz(i?). To extend 
this method to elements of degree > nz(i?), we proceed as follows. Let a ■ b be a 
monomial in S and suppose that \b\ = nz(i?). We use the normal form algorithm 
for low degree elements to obtain Nf(fe, B) = b\e\ + • • • + b m e m . To determine 
Nf (ab, B), write a = x^ 1 ■ ■ ■ x„ n . It is straightforward to check that 

m(ab,B) = (b 1 ,...,b m )A a 1 1 ■■■A a n "(l^e 1 ,...,l^e m ) t . 

Thus, the arithmetic complexity of the normal form algorithm is 0(|a|m 3 ) if 
one uses naive matrix multiplication or 0(|a|m 2 ' 376 ) if one uses state of the art 
methods [J]. To this one needs to add the complexity for computing Nf (b, B). 

Example 4.1. Suppose that we want to compute the normal form of x 17 with 
respect to the ideal / from Example 13.21 We have seen that {[x(y + z) 1 }, [y(y + 
z) 1 }, [z(y + z) 1 ]} forms a Ik-basis for R/I and that 

j II ON 

*--\\ \ -\) 

Since [x] is a basis element in degree one, the normal form of x 17 equals 
(1, 0, 0)A 16 (x(y + z) 16 , y{y + z) 16 , z(y + z) 16 )*. 
Since A 2 = 2 • A, wc have A 16 = 2 15 A. Hence 

Nf (a; 17 , x(y + z) w , y(y + z) 16 , z{y + z) w ) = 2 15 x{y + z) 16 . 
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If we know the variety of /, then the normal form computation can be 
simplified, see Example 14.41 



4.2 Upper bound of the elements in a Grobner basis 

To give an upper bound of the maximal degree of an element in a Grobner basis 
with respect to an ideal of projective dimension zero, we will use Gotzmann's 
persistence theorem. 

Recall that if h and i are positive integers, then h can be uniquely written 
as a sum 

'"'M:-D + - + (f 



h 

where 



rii > n.j_i > • • • > rij > j ' > 1. 



See |18j for an easy proof. This sum is called the binomial expansion of h in 
base i. Define 



h <l> = 



rii + 1\ f ri,_i + 1\ /rij + 1 



Before stating Gotzmann's theorem, recall that the Hilbert function of a 
graded algebra R is the map d i— ► dimj^i?^). 

Theorem 4.1 (Gotzmann's persistence theorem j§ ]). Let Hf be the Hilbert, func- 
tion of k[x±, . . . , x n ]/I, for any homogeneous ideal I. Let t denote the maximal 
degree of the generators of I . Then Hf(c?+1) = Hf(e?) <d> for some d> t implies 
thatm(d + 2) = Hf(d + l) <d+1> ,Hf(d + 3) = Hf (d + 2) <d+2> and so on. 

In the case of projective points, we have Hf(e? + 1) = Hf(e?) = m when 
d > m(R), thus we have Hf(d + l)< d + 1 > = Hf(d) <d> = m for d > nz(R) and 
hence 

N<d> _ f d + A , , fd-(m-2] 



W(d) 

\d+lj \d-(m-2) 

A Lex-segment set Ld on {x%, . . . , x n } is the \Ld\ biggest monomials of de- 
gree d in k[xi, . . . ,x n ] with respect to the lexicographical ordering. When L 
is a collection of Lex-segment sets, let L{L) denote the ideal generated by the 
elements in the Lex-segment sets. We call L(L) a Lex-segment ideal. When I is 
a homogeneous ideal, let | in(/)^| denote the number of monomials outside in(J) 
of degree d. Notice that | in(J)^| is independent of monomial ordering. 

Let / be a homogeneous ideal generated in degree less than or equal to d and 
let I be a collection of Lex-segment sets with maximal degree d. A property 
among Lex-segment ideals is that they have minimal growth (or maximal co- 
growth), in the sense that if | in(/) c | and | in(/(L)) c | agrees until degree d, then 

| in(I(L)) c d ,\ > | m{L< d ) d ,\ for all d' > d. 

See for instance |18j . 

Theorem 4.2. Let I be an ideal of projective dimension zero. A bound for the 
maximal degree of an element in a reduced Grobner basis is max(nz(_R), m). 



10 



Proof. Let d = max(nz(i?), to). Suppose that L is a collection of Lex-segment 
sets of degrees less than or equal to d, such that | in(/(L)) c | agrees with | in(/) c | 
until degree d. We then have 

m=\m(I(L)y d ,\ > \'m(I<d)d>\ > I Hhd')%\ = ™. 

This implies that | in^/^)^] = |in(7<d')^, | for all d 1 > d and hence there can 
not be any Grobner basis element of degree greater than d. □ 

This theorem is a generalization of the result in [1 ] , where it is shown that 
the last degree element of a Grobner basis is to in the case when I is unmixed. 
The bound in Theorcm l4.2l is sharp. Indeed, in Example l3.2l nz(i?) = 1, to = 3 
and a reduced Grobner basis with respect to DegRevLex had a generator in 
degree three, while in Example 14.21 below, we will see that nz(i?) = 3, m = 1 
and a reduced Grobner basis with respect to DegRevLex is {xy— z 2 , x 2 — xz, y 2 — 
z 2 , xz 2 — yz 2 , —yz 2 + z 3 }. 

4.3 Computing (R, A, I) given the ideal 

Suppose that we are given an ideal by its generators and that we know that 
diruis (Rd) = dimi 5 (i? £ i+i) for some d. What conclusions can be made from this in- 
formation? Unfortunately, not many. We do not know the dimension — indeed 
— the rings k[x, y, z]/(xy, yz, xz),k[x, y, z\j (x 2 , y 2 , z 2 ) and k[x, y, z]/(x 2 ,xy, xz) 
all have Ik-dimension three in degrees one and two. The first ring is of projective 
dimension zero and postulates in degree one. The second ring is artinian, while 
the third ring is of projective dimension one. However, we have the following 
simple observation. 

Lemma 4.3. Let I be a graded ideal in S and suppose that there is an element 
[/] G Ri such that [f]Rd = Rd+i- Then R is either artinian or of projective 
dimension zero. 

Proof. The ring S/(I + (/)) is artinian, hence S/I is of at most projective 
dimension zero. □ 

Lemma 4.4. Suppose that (/i,..-,/ n ) — I is generated by n elements in 
h[xo, . . . , x n ] and that there is an element [/] € Ri such that [f]Rd = Rd+i- 
Then R is of projective dimension zero. 

Proof. The ring S/(I + (/)) is artinian, hence (/i, . . . , /„, /) forms a regular 
sequence. But also (f\, ...,/„) forms a regular sequence, so S/I is of projective 
dimension zero. □ 

Even if we know that R is of projective dimension zero, it is also hard to 
tell whether or not the maximal ideal is associated. The following example 
shows that although dimk(-Rd) = dim^i^-i-i) and there is an element I such 
that [l]Rd — Rd+i, it does not hold that d > nz(R). 

Example 4.2. Let I = (x 2 -xz, xy- z 2 , y 2 ~z 2 ). Then Es(R, t) = 1 + 3t + 3t 2 + 
t 3 + t A + ■ ■ ■ and nz(i?) = 3. We have I = ( x — y, x — z) Q (z 2 ,y 2 , xy, x 2 — xz), 
V(I) = V((x — y, x — z)) = (1 : 1 : 1) and \J (z 2 , y 2 , xy, x 2 — xz) — (x, y, z). We 
can choose {[x], [y], [z]} och {[xz], [yz], [z 2 ]} as k-bases in degree one and two 
respectively and thus, the map from Ri to R2 induced by multiplication by [z] 
is injective. 
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Fortunately, as the next theorem shows, if we are only interested in comput- 
ing the variety, it is enough to find an I such that [l]Rd — Ra+i- 

Theorem 4.5. Let I be any homogeneous ideal and let R — S/I. Suppose that 
there exists an element [I] such that [l]Rd — Rd+i- Let [ft], ■ ■ ■ , [ft] be a k-basis 
for Rd+i- Let [ei], . . . , [et] be such that [e,-] [I] = [fi] . Let Aq, ■ ■ . , A n be such that 
Ai corresponds to multiplication with xi with respect to the bases [ei], . . . , [e*] 
and [fx], ...,[/*]. 

Suppose that p € V(I). Then e(p)* — (ei(p), . . . , et(p))' is a common eigen- 
vector to the Ai 's. Let \ be the eigenvalue of Ai corresponding to e(p) . Then 
p = (A : Ai : • • • : A„). 

Proof. By the definition of the matrix Aj we have 



[xj][e k ] = a$[lei] + ---+a$[le t ]. 

Thus, we get 

Xj(p)e k (p) = a$l(p) ■ ei{p) + ■■■ + a { k \ } l(p) ■ e t (p), 
or put in matrix form 

x^eipY = Ajl^eip) 1 . (1) 

Now e(p) can not be the zero vector, since otherwise we would have p € V(I + 
(ei)+- • - + (et)), which is a contradiction since S ' j (I + (ex) + • ■ - + (et)) i s artinian. 
With the same argument, l{p) must be non-zero. 

Hence ^j^y is an eigenvector of Aj with the eigenvalue . The theorem 
follows since 

(A : Ar : • • • : \ h ) = (x Q (p)/l(p) : n(p)//(p) : ■ • ■ : x n (p)/l(p)) 
= (xo(p) ■ xi(p) : ■■■ : x n {p)). 

□ 

Example 4.3. Let L = (y 2 , z 2 ,xz,xy). Then dimk(_Ri) = 3 and dimk(i?2) = 
2. We can choose [x], [y], [z] and [x 2 ], [yz] as k-bases in degrees one and two 
respectively. It is clear that [x + z\Ri = R3 and thus, by Lemma FOl we know 
that / is of at most projective dimension zero. We have [x + z][x] = [x 2 ] and 
[x + z] [y] = [yz] , so with respect to the bases [x] , [y] and [x + z] [x] , [x + z] [y] , we 
get 

Ax = (J 0) ' Ay = (0 0) and Az = (0 1 

There are two common eigenvectors for these matrices — (1>0) and (0, 1). The 
associated eigenvalues are 1,0,0 and 0,0,1 respectively. By Theorem 14. 5 [ we 
know that V{L) C {(1 : : 0), (0 : : 1)}. We have y 2 ((l : : 0)) = z 2 {{1 : 
0:0)) = xz((l : : 0)) = xy((l : : 0)) = 0, but z 2 ((0 : : 1)) ^ 0, so 
V(I) = {(1:0: 0)}. Thus, the second point was "false". 



Theorem 14.51 could also be used to compute the variety in Example l4.2l We 
leave this computation as an exercise to the reader. 

So suppose that we want to compute the variety of an ideal / which we 
suspect is of projective dimension zero. We propose the following procedure for 
a field k with enough elements. 
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Kl Compute the Grobner basis elements of degree 1, 2 and so on until we reach 
a degree d such that | in(/)^| > in(/)^ +1 | (this is the same as dim^-Rd) > 
dim k (i? d+ i) = t). 

K2 Choose a linear form [I] at random and check if Rd[l] = Rd+i- If it was 
not, choose another I. If we did not find such an element even after many 
tries, go back to stage Kl and compute more Grobner basis elements. 

K3 Choose a basis •••)[/*]} = [in(/)^ +1 ] f° r Rd+i and let e\, . . . ,e t be 

such that [eJ[Z] = 

K4 Compute the multiplication matrices with respect to {[et], . . . , [e*]} and 

{[/l], •••,[/*]}• 

K5 Determine a set of common eigenvectors for these matrices, either by using 
symbolic or numerical methods and use Theorem 14.51 to determine the 
variety of /. 

We refer the reader to the book |T5] and the citations therein for techniques 
to compute common eigenvectors using numerical methods. 

4.4 Computing (R, A, I) from the points 

Given a set of projective points P = {p\, . . . 7 Pm} one can form the vanishing 
ideal I(P), which consists of all polynomials vanishing on all of the points in 
P. The Hilbert series of R/I(P) is well studied but not completely understood, 
cf. [5]. The most common way of computing Hilbert series of an ideal defined 
by projective points has been studied by means of the projective Buchberger- 
Moller algorithm [TJ [14] . This algorithm computes a Grobner basis of a van- 
ishing ideal by computing a k-basis for the Ik-spaces Rq, R\ , . . . , until degree 
max(m, nz(i?)) and reducing potential Grobner basis elements with respect to 
this basis. We will present a reduced version of the projective Buchberger-Moller 
algorithm which instead of computing the Grobner basis of I(P) computes the 
triplet (R, A, I) and we will show that the behavior of our method is asymptot- 
ically better than the classical Buchberger-Moller-algorithm. 

Recall that we suppose that the representation of each projective point is 
fixed so that we can define evaluation of projective points in a unique way. 

A nice way to compute normal forms with respect to vanishing ideals of 
projective points is by evaluation: Given a form / of degree d and a vector 
space basis e\, . . . , e m of Rd, we obtain the normal form ct\e\ + • • • + a m e m , 
where the ai's are chosen to satisfy f(pi) — ceiei(pi) + ••• + a rn e m (pi) for 
i = 1, . . . , m. The normal form does not depend on the choice of representation 
of the points. Computing normal forms by means of evaluation is the key engine 
behind the graded Buchberger-Moller algorithm and the variation of the method 
given below. 

When studying ideals of projective points, one can always assume that n + 
1 < m. Indeed, we have the following lemma, which is a graded version of 
Lemma 5.2 in [T2"] . 

Lemma 4.6. Let E = {xi , . . . ,Xi w } be any subset of the variables such that 
E(P) and {xq, . . . , x n }(P) has same rank. Let n be the natural projection from 
P n (k) to P n (k) defined by 7r((ao : • • • : a n )) = (fflj : • • • : a;_). Then q±, . . . , q m 
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are distinct where qi — Tr(pi). Moreover, with Q = {q\, . . . , q rn } and with R — 
k.[xi , . . . , Xi—]/I{Q), the graded algebras R and R are isomorphic. 

Proof. Suppose that Xi £ E. Then Xi(P) = a\Xi 1 (P) + • • • + ctnXi-(P). Hence 
Xi — a\Xi t + • • • + otfiXi- 6 /. Since i?>i is generated in degree one, it is clear 
that the elements in E generates R>\. Since the evaluation on q^s and the p^'s 
agrees on the elements ] , it follows that R and R are isomorphic 

as graded algebras. It is the clear that q±, . . . , q m are distinct. □ 

Remark 2. In a more subtle way, Lemma \4-.6\ actually follows directly from the 
projective Buchberger-Moller algorithm. 

We now give a variant of the projective Buchberger-Moller algorithm for 
building the triplet (R,A,l) from the points. As for the Buchberger-Moller 
algorithm, this algorithm is based on the evaluation method to compute normal 
forms. But it differs from the Buchberger-Moller algorithm in the sense that 
it is focused on giving the multiplication tables with respect to the variables 
rather than giving a Grobner basis for the ideal. 

LI Compute a non-zero divisor / of degree one by using the method in Propo- 
sition 13.21 

L2 Initiate the lists P>o = Lq = [1] and Initials = []. Let d = 0. 

L3 If rank(B d (P)) = \P\, then nz(R) = d. Return B ,...,B d and I. Other- 
wise, increase d by one, let Bd = [] and let Ld be the list of all monomials 
of degree d which are not multiples of an element of Initials. 

L4 If Ld is empty, go to step L3; otherwise choose the monomial t — min^ (L<i) 
with respect to a fixed monomial order and remove it from Ld- 

L5 If t(P) can be written as a linear combination of the rows in Bd(P), then 
add t to the set Initials and continue with step L4. Else, append t to Bd 
and continue with step L4. 

The correctness of the method is a direct consequence of the projective BM- 
algorithm, since the sets Bq, . . . ,Bd are computed in the same way using the 
two methods. Thus Initials will generate in(J)^, while B will be the complement 
of in(J)^. By using another selection method in step L4, it is possible to obtain 
a basis which is not necessarily the complement of an initial ideal (It is an easy 
exercise to check that the termination of the algorithm does not depend on the 
selection method). 

We implicitly assume that we have used Lemma 14.61 so that n < m. This 
preprocessing can be done using 0(nm 2 ) arithmetic operations, since we test for 
linear dependence n times. It is straightforward to lift the result in and show 
that the complexity of the algorithm is dominated by the arithmetic operations 
and not the monomial manipulations. The number of arithmetic operations for 
the step LI is bounded by 0(m 3 ) by an elementary analysis of the method in 
Proposition ^. 21 For each degree d during the algorithm, we need to check linear 
dependence at most min(m, n)m times. Thus, the arithmetic complexity of the 
method is a cost of at most 0(post(i?) min(m, n)m 3 ). The original analysis [14] 
of the Buchberger-Moller algorithm reports the complexity 0(nm 4 ). (Although 
it is possible to show that the performance is 0(min(m, n)nm 3 ).) 
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Example 4.4. In Example 3.5 in [J, the point set P = {(0 : 2 : 5), (0 : 1 : 2), 

(1 : 3 : 1), (4 : 3 : 4), (2 : 5 : 4), (1 : 4 : 4)} is considered. The Grobner basis 
with respect to DegRevLex and x -< y -< z is generated in degrees three and 
four and Rs{R, t) = 1 + it + 6t 2 + 6t 3 + ■ ■ ■ . 

With our approach, we would first fix the coordinates P — {(0, 2, 5), (0, 1,2), 
(1,3,1), (4,3,4), (2,5,4), (1,4,4)} and then compute L = {l},B a = {1} and 
L\ = {x,y, z}. Since {x, y, z}(P) has full rang, we will have B\ = {x,y, z}. 
Thus, Initials is empty and L 2 = {z 2 ,yz,xz,y 2 ,xy,x 2 }. It turns out that also 
{z 2 , yz, xz, y 2 , xy, x 2 }(P) has full rang, so B2 = {z 2 ,yz,xz,y 2 ,xy,x 2 }. Since 
\P\ = 6, the algorithm stops and we know that Hs(R,t) = 1 + 3t + 6t 2 + 
6t 3 + • • • . It is immediate that y{pt) 7^ for i — 1, ... ,6, so we can use 
-B3 = {yz 2 , y 2 z, xyz, y 3 , xy 2 , x 2 y} as a k-basis in degree three, and in general 

B d = {y d - 2 z 2 , y^z, xy^ 2 z, y d , xy d -\x 2 y d - 2 }. 

Say that we want to compute the normal form of x e + z 6 . If we do it by 
evaluation, we solve the linear equations 

(a; 6 + z 6 )(pi) = (aiy 4 z 2 + a 2 y 5 z + a 3 xy 4 + a^y 6 + a 5 xy 5 + a e x 2 y 4 )(p l ) 

for i = 1, . . . , 6 which is equivalent to perform Gaussian elimination on a (6 x 6)- 
matrix. As result, we get 

^r, r fi n , 2083926583 R 11603225231 . 8111541583 = 

Nf(x + z B) = v xv v z 

v ' ' 23522400 y 470448000 y 26136000 y 

327280970021 2 4 17527852333 4 2 127511218609 4 
940896000 X V + 117612000 V Z + 313632000 Xy Z ' 
We could also use the multiplication matrices. Notice that A x ,A y and A z 
share the six linear independent eigenvectors 

(z 2 (p l ),yz(p l ), xz(pi), y 2 (p t ), xy(pi), x 2 {p l )) t . 

Thus, if we let T = B 2 {P) t , we have A x = TD X T~ X , A y = TDyT^ 1 and A z = 
TD Z T~ X , where D x = diag(0,0, 1/3, 4/3, 2/5, 1/4), D y = diag(l, 1, 1, 1, 1, 1) and 
Z> z = diag(5/2,2,l/3,4/3,4/5,l). 

So to compute the normal form of x 6 + z 6 , we can start by computing 
M(x 2 ,B) = x 2 = (0,0,0,0,0, 1) ■ (z 2 ,yz,xz,y 2 ,xy,x 2 ) and M(z 2 ,B) = z 2 = 
(1,0,0,0,0,0) • (z 2 ,yz,xz,y 2 ,xy,x 2 ). We get Nf(a; 6 + z 6 ,B) = Nf(x 6 ,S) + 
Nf^ 6 ,^), where 

Nf (x 6 , B) = (0, 0, 0, 0, 0, l)Ai(y 4 z 2 , y 5 z, xy 4 z, y 6 , xy 5 , x 2 y i ) t 

= (0, 0, 0, 0, 0, \)TD x T-\y 4 z 2 , y 5 z, xy 4 z, y 6 , xy 5 ,x 2 y i ) t 

and 

Nf (z 6 , B) = (1, 0, 0, 0, 0, 0)A 4 z (y 4 z 2 , y 5 z, xy 4 z, y 6 , xy\x 2 y 4 ) t 
= (1, 0, 0, 0, 0, 0)TZ£T- Vz 2 , y 5 z, xy 4 z, y 6 , xy 5 , x 2 y 4 f. 
Of course, we get the same result as above. 
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4.4.1 Computing separators 

A family of separators with respect to a set of affine points P = {pi, . . . ,p m }, is 
a set {/i, . . . , f m } of polynomial functions such that fi(pi) = 1 and fi(pj) = if 
i 7^ j. It is easy to see that the separators forms a k-basis for k[xi, . . . , x n ]/I(P). 

When the points are projective, we say that {/i, . . . , f m } is a set of separators 
if fi(Pi) and fi(pj) = when i ^ j. When all separators are of the same 
degree d, they constitute a k-basis for Rd- If k contains at least m elements so 
that there exists a non-zero divisor I, we can construct a separator-k-basis for 
Rd+i as fil 1 , . . . , f m l l . 

In [12] , two methods for computing separators with respect to a collection of 
affine points are discussed. It is possible to lift this method to the projective set- 
ting. Both methods perform the same number of arithmetic operations. We will 
illustrate one of the methods by an example. In Z 5 , consider p\ = (1, 2, 0, 1, 1), 
p 2 = (1, 0, 1, 1, 2), p 3 = (1, 2, 0, 3, 3), Pi = (0, 0, 2, 0, 4), p 5 = (0, 0, 2, 1, 5) and 
Pq = (2, 1,3, 1,6). We associate the following table to this point set 



/ 1 


1 


1 








2 


Si = 


{{1,2, 3}, {4, 5}, {6}} \ 


2 





2 








1 


s 2 = 


{{1,3}, {2}, {4, 5}, {6}} 





1 





2 


2 


3 


s 3 = 


{{1,3},{2},{4,5},{6}} 


1 


1 


3 





1 


1 


£4 = 


{{1}, {3}, {2}, {4}, {5}, {6}} 


V i 


2 


3 


4 


5 


6 


£5 = 


{{1}, {3}, {2}, {4}, {5}, {6}} J 



The sets on the right hand side are also described by an example: The set {1, 3} 
on the second row shows that p\ and p^ agree on the first two coordinates. 
When computing such a table from a point set, one obtains a matrix aj, where 
Cij is the first position where pi and pj differ. 

This matrix is used to compute the separators and it is clear that 



n 



Pica 



Pi, 



Pjc i: 



satisfies Qiijpj) = if i ^ j and Qi(pi) = 1. 

It is showed in |12j that at most nm + m 2 arithmetic comparisons are used 
to compute the matrix . (In fact a slightly improved upper bound is given.) 

We will now show how to make use of the matrix Cy to compute projective 



separators. If we let Sij{jpi) 7^ and Sij(pj) - 
projective separators for p\, . . . ,p m , where 

Qi = Sij . 



0, then Qi, . . . , Q m is a set of 



(2) 



Suppose that each point pi is normalized in the sense that the first non- 
zero position equals one. It is then clear that we can use the affine method to 
compute the matrix (cy) with respect to the points. 

We will now give an explicit formula for each Sij . To simplify notation, let 
h = CjA. 



• If pih = 0, then pjh 7^ 0. Let h! be the least position such that pw 
and let S^ = PjhXh' ~ Pjh'Xh- 

• Else, if p^ 7^ but pjh = 0, then let = Xh- 
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• Finally, suppose that p^ ^ and pjh ^ 0. Since Pi and pj agrees on 
all coordinates less than h and p^ ^ Pjh, there is a h' < h such that 
Pih' = Pjh' = 1- Thus, let Sij = PjhXh' - Pjh'Xh = P hXh> - x h . 

Notice that we can choose the index h! occurring in the two situations as the 
first entry where pi equals one. It is clear that we can determine the first non- 
zero index of each point using at most nm arithmetic comparisons. We have 
proved the following theorem. 

Theorem 4.7. Let P = {pi, . . . ,p m } be a set of distinct projective points. Then 
we can compute a set of separators of degree m — 1 with respect to P using at 
most nm + m 2 arithmetic operations. 

Example 4.5. Let Pl = (1 : 2 : : 1 : 1 : : 3 : 5), p 2 = (1 : : 1 : 1 : 2 : : 3 : 
5), p 3 = (1 : 2 : : 3 : 3 : 1 : 2 : 0) and let Pi = (0 : 1 : 1 : : 2 : : 1 : 0). We 
will show how to compute Q\. We have ci 2 = 2 and p\2 = 2 and P22 = 0. Thus, 
S12 = X2- We have C13 = 4 and pu = 1 and P34 = 3. Since pu = p 3 i = 1, we 
let S13 = P34X1 — X4 = 3xi — X4. We have C14 = 1 and p 4i = 0, so Si 4 = x\. 
Hence Qi = x 2 (3xi — x 4 )xi. 

5 Discussion 

In a forthcoming paper we will generalize parts of the results to rings of arbitrary 
projective dimension. 
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